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Abstract 

O I We construct the equation of state (EOS) of nuclear matter at finite temperature and 

^ ' density with various proton fractions within the relativistic mean field (RMF) theory for 

I'-j • the use in the supernova simulations. We consider nucleus, alpha-particle, proton and 

^ , neutron in equilibrium at densities smaller than about p_B ~ 10^^'^g/cm'^ by minimizing 

the free energy of the nuclear matter. The calculation is based on the RMF theory with 
the parameter set, TMl, which has been demonstrated to provide good accounts of the 
ground and the excited state properties of finite nuclei. We tabulate the outcome for 
various densities as the pressure, the free energy, entropy etc, at enough mesh points in 
the density range ps = 10^'^ ~ 10^^'^ g/cm^ and the temperature range T = ~ lOOMeV 
and the proton fraction range of 1^ = ~ 0.56 to be used for the supernova simulations. 
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1 Introduction 



The fate of massive stars with mass of more than 8 times the solar mass is spectacular in 
the end when they have used up their nuclear energy by creating the Fe core, which has the 
largest binding energy per particle. When the Fe core mass exceeds the Chandrasekhar mass 
supported by the degenerate electron pressure against the gravitational force, the Fe core starts 
to collapse until the central density becomes around that of the normal matter density. There 
the strong repulsive force and the Pauli effect among nucleons work against the gravitational 
collapse and swing back the matter motion outwards. If this outward motion is strong enough 
to blow off the outer layers, the stellar explosion is observed as supernova. This is called the 
prompt explosion. The strength of the outward shock motion depends crucially on the equation 
of state (EOS) of nuclear matter. 

The general belief in the recent years is, however, that the outward shock propagation 
alone is not strong enough to succeed supernova explosion. Most of the gravitational collapse 
energy is carried away by neutrinos. If some small fraction of the neutrino energy is dropped 
in the stellar matter just behind the slowing down shock region with possible help of the 
hydrodynamic instability, the outward shock motion can survive to blow off the outer materials. 
This is called the delayed explosion. Whichever the mechanism of the supernova explosion is, 
we need detailed knowledges on the EOS of nuclear matter and also the neutrino reaction rates 
in nuclear matter in the wide ranges of density and temperature. We have started a project 
to provide these informations, the EOS and the neutrino reaction rates, by studying carefully 
the existing knowledge on finite nuclei including the unstable nuclei and the weak processes. 
In this paper, we discuss the results on the EOS. 

The knowledge on nuclear matter necessary for supernova simulations is much more than 
that obtained from the study of finite nuclei. The density may go up to 5 to 10 times the 
normal matter density in the stellar phenomena. The temperature may climb up to 30 to 50 
MeV. The proton fraction may go down to zero, the pure neutron matter. Hence, we have to 
make a large extrapolation to all directions from our knowledge on nuclear physics which are 
available from laboratory experiments. Hence, we have to construct a theoretical model, which 
has a strong basis from the microscopic theory and at the same time is easy to handle for many 
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observables. This theoretical model has to be checked with many experimental facts on finite 
nuclei. 

In recent years, there was a very important achievement on the theoretical understanding 
on nuclear physics. The use of the relativistic many body theory enabled us to remove the 
long standing problem on the nuclear saturation property, which is not reproduced in the non- 
relativistic approach . The relativistic theory can reproduce many observables on the nuclear 
structure including the unstable nuclei and also the nuclear reactions, in particular the spin 
observables 0, 0, ^, ^. The relativistic Brueckner-Hartree-Fock (RBHF) theory reproduced 
the saturation property of the nuclear matter very close to the 'experimental' value and re- 
moved the disease of the non-relativistic Brueckner-Hartree-Fock theory, which provided the 
saturation points lying on the Coester line. The relativistic framework contains very strong 
density dependent repulsive component due to Z-graph three body process in the language of 
non-relativistic framework. The strong spin-orbit force and the spin dynamics are the natural 
consequences of the relativistic equation of motion. 

Hence, we chose the relativistic mean field (RMF) theory, which is relatively easy to use and 
at the same time carries the essential features of the microscopic theory as the RBHF theory, 
as a phenomenological model to describe nuclei and nuclear matter at various circumstances 
ip. First, the RMF lagrangian was constructed to provide the results of the RBHF, for which 
we include the non-linear uo meson coupling term in addition to the generally used non-linear 
(J meson coupling terms. We fixed the model parameters as the coupling constants and the 
masses so as to reproduce many key experimental data on the ground state masses and the radii 
of nuclei including unstable ones 0]. The parameter set was named TMl and then applied to 
many observables as deformed nuclei, triaxial deformations, nuclear excitations |P, 0, §. The 
agreement with the existing data is extremely satisfactory. Hence, the RMF theory with the 
TMl parameter set could go through the critical check from Nuclear Physics. The RBHF 
results for the symmetric nuclear matter were also compared very nicely up to high density. 

This was the place we decided to go on to construct the EOS table for the astrophysical 
use. We want to state here our philosophy. The purpose of making this table with only the 
nuclear degree of freedom is to work out the best EOS with the known experimental and the 
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theoretical knowledge on nuclear physics. This EOS should be considered as the basis EOS. 
At high density and temperature, however, there may be the need of the mixture of hyperons, 
various mesons and even quarks. As our knowledges on these degrees of freedom become more 
solid, we should include them for the EOS table, in particular at high density. 

We work out the EOS of nuclear matter in the wide parameter ranges in density, temperature 
and proton fraction. The density range is Pb = 10^'^ ~ IQ^^ '^g / cm? , the temperature range 
is T = ~ IQQMeV and the proton fraction is 1^ = ~ 0.56. In very high density and 
high temperature region, the nuclear matter is homogeneous and hence the calculation is not 
difficult. However, when the density is below ps ~ l/3po the matter becomes inhomogeneous. 
In order to treat the inhomogeneous matter, we take the Thomas-Fermi approximation, where 
the inhomogeneous matter is composed of a lattice of spherical nuclei with protons, neutrons 
and alpha-particles 0. We then minimize the free energy of the nuclear system with respect to 
the proton, neutron and alpha-particle densities in the nuclear interior and the exterior regions. 
This procedure has to be performed for each temperature, density and the proton fraction, which 
is extremely cumbersome. At the same time, we have to repeat the minimization procedures 
with good accuracy in order to get the smooth physical quantities as the free energy, pressure 
and the entropy for the use of supernova simulations. Some part of the present research project 
is published in Ref.[0. The published materials do not contain the contributions of the alpha- 
particles. The new part of this paper is the inclusion of the alpha-particles and hence we shall 
make detailed comparisons of the outcome with respect to the alpha-particle contributions. 

This paper is arranged as follows. We describe the RMF theory and the TMl parameter 
set in Sect. 2 and the formalism of the Thomas- Fermi approximation with the alpha-particles 
in Sect. 3. The numerical results are discussed in Sect. 4. The conclusion of this paper is stated 
in Sect. 5. We describe the definitions of the physical quantities tabulated in the EOS table in 
Appendix A. 
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2 Relativistic mean field theory 

We adopt the RMF theory with non-hnear a and uj terms to describe homogeneous nuclear 
matter [Q. We start with the lagrangian given by 

JC^rmf = i'iilt.d^' - M - g^a - g^-i^uj^ - Qp-i^Tap''^']^ (1) 

1 111 

H — dncrd^a mla'^ oacr^ q^a'^ 

2 M 2 3 4 

Here, if) denotes a SU(2) baryon field of mass M (proton and neutron); cr, cu'^, and p"'^ are 
0", cij, and p meson fields with masses mcr,mi^, and rrip, respectively. W^'^ and R"-'^'^ are the 
antisymmetric field tensors for tu^ and p"^, which can be written as 

W^"" = d^'uj^ -d'^uj^, (2) 
^a/.^ = SV"" - + ^pe"VV- (3) 

In the lagrangian, the constants gaiQui-, Qp are the coupling constants for the interactions between 
mesons and nucleons, the coefficients g2 and g^ are the self-coupling constants for a meson field, 
while C3 is the self-coupling constant for uj meson field. It is known that the inclusion of the non- 
linear 0" terms is essential to reproduce the properties of nuclei quantitatively and a reasonable 
value for the incompressibility, while the non-linear uj term is added in order to reproduce 
the density dependence of the vector part of the self-energy of the nucleon obtained in the 
RBHF theory. The lagrangian contains the meson masses, the coupling constants, and the self- 
coupling constants as parameters. We adopt the parameter set TMl listed in Table 1, which 
was determined in Ref . as the best one to reproduce the properties of finite nuclei in the wide 
mass range in the periodic table including neutron-rich nuclei. The RMF theory with the TMl 
parameter set was also shown to reproduce satisfactory agreement with experimental data in 
the studies of the nuclei with deformed configuration and the giant resonances within the 
RPA formalism With the TMl parameter set, the symmetry energy is 36.9 MeV and the 
incompressibility is 281 MeV. 
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Starting with the lagrangian, we derive a set of the Euler-Lagrange equations. The Dirac 
equation for the nucleon field is given by 

(^7^9^ -M-g„a- g^^^uj^ - g.^^.Tap'"') = 0, (4) 

and the Klein-Gordon equations for the meson fields are given by 

d^d'^a + mla = -ga'ipi) - g2a'^ - g3(^^, (5) 
d^W'^ + mluj'' = g^i,ri^-c:,{u,uj'')uj^, (6) 
d^R'^^ + my^ = gpt/jTar^ + gpeabcPlR'"^. (7) 

These equations are coupled nonlinear quantum field equations, which are very difficult to solve 
exactly. We adopt the relativistic mean field approximation as described in Ref. |TT|, in which 
the meson fields are treated as classical fields, and the field operators a, uj'^, and p"^ are replaced 
by their expectation values (cr), (uj^), and {p"-^). In the present study, we consider first static 
infinite matter so that we obtain simplified equations, where the derivative terms in the Klein- 
Gordon equations vanish automatically due to the translational invariance of infinite matter. 
The spatial components of the vector meson fields vanish under the rotational symmetry. For 
the isovector- vector meson field p"^, only the third isospin component has a nonvanishing value 
because of the charge conservation. Hence, the equations for meson fields are reduced to 

= i^) = -^^{i'^)-^^{92^o+93a',), (8) 

coo = (c.°) = 4(^7V)-^C3o;o^ (9) 
mi mi 

Po ^ (p'°) = 4(V^^37V). (10) 
ml 

The Dirac equation then becomes 

(-iafcV'' + (3M* + g^ujQ + gpT^po) ipis = Sisipis, (11) 

where the index i denotes the isospin degree of freedom (proton and neutron) and s denotes the 
index of eigenstates of nucleon, M* = M + g^a^ is the effective mass, Sis is the single-particle 
energy. 

Nucleons occupy single-particle orbits with the occupation probability fis. At zero temper- 
ature, fis = 1 under Fermi surface, while fis = above Fermi surface. For finite temperature. 



the occupation probability is given by Fermi-Dirac distribution, 

1 1 



fii 



ft! 



l + exp[{eis- /T] i + exp 
1 



P + -^i] /T 
1 



l + exp[{-eis + fJ,i) /T] 



exp 



f.2 + + Ui] /T 



(12) 
(13) 



Here the index i and i denote nucleons and antinucleons, and are the single-particle 
energy of nucleons and antinucleons, respectively. The relation between the chemical potential 
fii and the kinetic part of the chemical potential Ui is give by 



The chemical potential /Xj is related to the number density of nucleon rii as 

ry poo 



(14) 



Hi 



dkk ifik-fik) 



(15) 



27r2 Jo 

where 7 is the degeneracy factor for the spin degree of freedom (7 = 2), the quantum number 
s is replaced by the momentum k when we do the integration in the momentum space instead 
of the summation of eigenstates. The equations of meson fields can be written as 



Po 



mfr" 27r^ Jo 



'o- i 



M* 



k^ + M* 



ifik + fik) 2 (^'2^0 + 93(^0 1 ' 



— [Hp + n„) 



— ^ [Up - Un) . 



m 



(16) 

(17) 
(18) 



Here, Up and n„ are the proton number density and neutron number density as defined in 
Eq.(15), and we denote = + n„ as the baryon number density. We solve these equations 
self-consistently. The thermodynamical quantities are described in Ref. [^, and we just write 
the expressions here. The energy density of nuclear matter is given by 



^ roc I 1 11 

E ^ X dk k' ^k^ + M*2 (/,, + /,,) + -mlal + -g^al + -g,a, 



(19) 



+9uj^^Q {np + rin) - ^mlujQ - ^CsUq 



+9pPo {np -Un)- -m^po, 
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the pressure of nuclear matter is given by 

ry POO -j^ -j^ -j^ 

^ 67T^Jo Jk^ + M*^ 2 3 4 



and the entropy density is calculated by 



s = E ^ / ^^^'^ i~fik In fik - (1 - fik) In (1 - fik) - fik In fik - (1 - fik) In (1 - fik)]- 

(21) 



3 Thomas- Fermi approximation 

In the range, T < 15MeV and pb < lO^^'^fyf/cm^, where heavy nuclei may be formed in order 
to lower the free energy, we perform the Thomas-Fermi calculation based on the work done 
by Oyamatsu 0. In this case, the non- uniform matter can be modeled as a mixture of free 
neutrons, free protons, alpha-particles, and a single species of heavy nuclei, while the leptons 
can be treated as uniform non-interacting particles separately. For the system with fixed proton 
fraction, the leptons play no role in the minimization of the free energy. Hence we mainly pay 
attention to baryon contribution in this calculation. Hereafter, we will not mention the leptons 
frequently, while we should keep in mind that there exists an uniform lepton gas existing 
everywhere. 

We assume that each heavy spherical nucleus is located in the center of a charge-neutral 
cell consisting of a vapor of neutrons, protons and alpha-particles. The nuclei form the body- 
centered-cubic (BCC) lattice to minimize the Coulomb lattice energy. It is useful to introduce 
the Wigner-Seitz cell to simplify the energy of the unit cell. The Wigner-Seitz cell is a sphere 
whose volume is the same as the unit cell in the BCC lattice. We define the lattice constant a 
as the cubic root of the cell volume. Then, we have 

= a^ = Nb/ub, (22) 

where Nb and ub are the baryon number per cell and the average baryon number density, 
respectively. We calculate the Coulomb energy using this Wigner-Seitz approximation and add 
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a correction energy for the BCC lattice 0. This correction energy is neghgible unless the 
nuclear size is comparable to the cell size. 

We assume the nucleon distribution functions nj(r) {i = n for neutron, i = p for proton) in 
the Wigner-Seitz cell as 



rii (r) 



1 



+ nr*, 0<r <R 



(23) 



Ri<r< Rceii 



where r represents the distance from the center of the nucleus and Rceii is the radius of the 
Wigner-Seitz cell defined by the relation, 

The density parameters n*" and n°^* are the densities at r = and r > Ri, respectively. The 
parameters Ri and ti determine the boundary and the relative surface thickness of the heavy 
nucleus. 

For the distribution function of alpha-particles na{r), which should decrease as r approaches 
the center of the heavy nucleus, we assume 

1 3 

I -nr^ 1 - ( -^V 

Ua (r) 



+ 0<r <R, 



, (25) 



R„<r <R, 



■p ^ I ^ J^cell 



which will give n^ir = 0) = and na{r > Rp) = Here we use the same parameters Rp 

and tp for both proton distribution function and alpha-particle distribution function in order 
to avoid too many parameters in the minimization procedure. The parameters i?„ and t„ may 
be a little different from Rp and tp due to the additional neutrons forming a neutron skin in the 
surface region. In order to illustrate these parameters, we show the distributions of nucleons 
and alpha-particles in Fig.l for the case at T = lOMeV, Yp = 0.3, and ps = 10^^'^ g/cm^. Here 
we define the baryon mass density as = rriunB with being the atomic mass unit. For 
the system with fixed temperature T, proton fraction Yp, and baryon mass density pb, there 
are eight independent parameters in the ten variables; a,n'^,n°^\ Rn,tn,n^j^,n°^^, Rp,tp,n'^^. 
The thermodynamically favorable state is the one that minimizes the free energy density with 
respect to these eight independent parameters. 
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In this model the free energy density contributed by baryons is given as 

/ = /a'= { - T SceU ) / a^ (26) 
where the free energy per cell Fceii can be written as 

Fceii = {Ebulk + Es + Ec) — TSceii = Fbulk + -E's + -E'c- (27) 
The bulk energy Ebuik^ entropy Scdu and bulk free energy Fbuik are calculated by 

Ebulk = / e{nn{r),np (r) , Ua (r) ) cfr, (28) 



cell 

Sccii = / s{nn{r) ,np (r) , (r) ) d^r, (29) 

J cell 

Fbulk = / firin (r) , Up (r) , (r) ) d^r. (30) 

Jcell 

Here e ( ra„ (r) , Up (r) , n^, (r) ), s ( n„ (r) , (r) , ra„ (r) ), and / ( n„ (r) , Up (r) , n^, (r) ) are the 
local energy density, entropy density, and free energy density at the radius r, where it can be 
considered as a mixed uniform matter of neutrons, protons, and alpha-particles. The local free 
energy density is given as the sum of the nucleon and alpha-particle contributions; 

/ ( Tin (r) , Up (r) , ria (r) ) = /at + /„, (31) 

where and fa are the contributions from nucleons and alpha-particles, respectively. We use 
the RMF theory to calculate the free energies of neutrons and protons when density is not 
extremely low (?7,„ + Up > 10^^ fm^^), while the classical ideal gas approximation is adopted 
at low density (n„ + np < 10^^ fm'^). The connection between the RMF results and the ideal 
gas approximation is very smooth. 

The alpha-particles are treated as non-interacting Boltzman particles, which is considered 
as a reasonable approximation because the number density of alpha-particles is not so large 
even though the average baryon density is high. The free energy density of alpha-particles is 
given by 

fa = -Tna[ln{8nQ/na) + l]-naBa, (32) 
where we have used the abbreviation uq = ^ , Ba = 28.3MeV is the binding energy of 



an alpha-particle relative to free neutrons taken from Ref.[12|. We have to take into account the 
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volume of alpha-particle, otherwise the alpha-particle fraction will become some big number at 
high densities, where the alpha-particles should actually disappear. After we take into account 
the volume excluded by the alpha-particles, the free energy densities of nucleons and alpha- 
particles become, 

/Ar(n„,np) = (1 - M)/Ar(n„,np), (33) 

fa{na) = {1-U)fa{ha), (34) 

where u = UaVa is the fraction of space occupied by alpha-particles, the effective volume of an 
alpha-particle, Va = 24/m~^, is taken from Ref.[^. We denote the effective number density of 
neutrons {i = n), protons {i = p), or alpha-particles {i = a) as fij = nj/(l — m). The inclusion 
of the volume excluded by the alpha-particles has negligible effect in the low density region, 
while it is necessary for the calculation at high density. 

As for the surface energy term Eg due to the inhomogeneity of nucleon distribution, we take 
a simple form as, 

Es= f Fo I V ( n„ (r) + (r) ) ^ d^r. (35) 

Jcell 

The parameter Fq = 70 MeV ■ fm^ is determined by doing the Thomas- Fermi calculations of 
finite nuclei so as to reproduce the gross properties of nuclear mass, charge radii and the beta 
stability hne as described in the appendix in Ref. [0. 

The Coulomb energy per cell Ec is calculated using the Wigner-Seitz approximation and 
add a correction term for the BCC lattice |p 

Ec = l f e K (r) + (r) - n,] <p{r)d\ + A^c, (36) 

2, J cell 

where 0(r) stands for the electrostatic potential calculated in the Wigner-Seitz approximation. 
He is the electron number density of uniform electron gas {ue = YpUB), AEc is the correction 
term for the BCC lattice, which can be approximated as 

AEc = Cscc^-^^, (37) 
a 

here a is the lattice constant as defined in Eq.(22), the coefficient Cbcc = 0.006562 is taken 
from Ref.[^, Znon is the non-uniform part of the charge number per cell given by 

3 

Anr^dr. (38) 
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Because of the long-range nature of the Coulomb interaction, the Coulomb energy will depend 
on the lattice type. This dependence was discussed in more details in Ref.[^. The system 
prefers the BCC lattice because the BCC lattice gives the lowest Coulomb energy. 

4 Equation of state of nuclear matter 

We construct the EOS table of nuclear matter in the wide ranges of the temperature T, proton 
fraction Yp, and baryon mass density as listed in Table 2 for the use of supernova simulations, 
and we also add the results for zero temperature case (T = 0) and pure neutron matter case 
{Yp = 0). The physical quantities we provide in the table and their definitions are described 
in Appendix A. The leptons can be considered as uniform non-interacting particles, which 
are relatively easy to deal with, so we provide the baryon part of EOS without the lepton 
contributions in this table. 

For each T, Yp, and pb, we have to determine the thermodynamically favorable state of 
nuclear matter, which has the lowest free energy in this model. We perform the minimization 
of the free energy for both non-uniform matter and uniform matter. Here the phase of heavy 
nuclei formed together with free nucleons and alpha-particles is referred to as non-uniform 
matter while the phase of nucleons mixed with alpha-particles without heavy nuclei is referred 
to as uniform matter. For non-uniform matter, the minimization procedure is realized by using 
the Thomas-Fermi method, which includes eight independent parameters as described in Sect. 3. 
For uniform matter, we do the minimization with respect to converting two protons and two 
neutrons into an alpha-particle, so there is only one independent parameter in this minimization 
procedure. By comparing the free energies of the non-uniform matter and uniform matter, we 
determine the most favorable state of nuclear matter at this T, Yp, and pb- 

We discuss first the phase transition between uniform matter and non-uniform matter. The 
phase transition at zero temperature has been discussed in Ref . ||T0| , which does not change 
after including alpha-particle, because both the proton fraction and the alpha-particle fraction 
at zero temperature are equal to zero. In Fig.2, we show the phase boundary dependent on the 
temperature T at 1^ = 0.01,0.1,0.2, and 0.4. The solid curve is the boundary of non-uniform 
matter phase where the heave nuclei are formed so that the heavy nucleus fraction is more 
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than zero denoted by Xa > 0. The dashed curve is the boundary where the alpha-particle 
fraction changes between < 10~^ and X^ > 10~^. The definitions of Xa and X^ are given 
in Appendix A. At low density, the thermodynamically favorable state is homogeneous nucleon 
gas with a small fraction of alpha-particles, the heavy nuclei are formed at some medium 
densities where the system can lower the free energy by forming heavy nuclei, then it becomes 
uniform matter as density increases beyond pb ~ 10^^'"^ g / cm^ . It is seen in this figure that the 
density where the non-uniform matter phase appears depends on the temperature very strongly, 
while the density where the non-uniform matter phase disappears is nearly independent of the 
temperature. As temperature increases, the density range of the non-uniform matter phase 
becomes narrower. The non-uniform matter phase disappears completely at high temperature. 
In Fig. 3, we show the phase boundary dependent on the proton fraction 1^ at T = 0, 1, 5, 
and lOMeV, this dependence is relatively weak except at very low Yp limit. The non-uniform 
matter phase disappears at lower temperature with a small Yp because it is difficult to form 
heavy nuclei with a small proton number. 

We plot in Fig. 4 the fractions of neutrons, protons, alpha-particles, and heavy nuclei as 
functions of the baryon mass density at 1^ = 0.2 with various temperatures. The alpha- 
particle fraction Xa increases as baryon mass density ps increases, but the formation of heavy 
nuclei brings rapid drops for the alpha-particle fraction and the nucleon fractions, because the 
heavy nuclei use up most of the nucleons. From this figure, we can find the alpha-particle 
fraction becomes smaller as the temperature increases. 

In Figures 5a and 5b, we display the nuclear mass number A and the proton number Z as 
defined in Appendix A as a function of the baryon mass density pb with various proton fraction 
1^ at T = IMeV. Both A and Z increase sharply before the phase transition. The nucleus 
becomes larger before the system turns into the uniform matter phase. This behavior is similar 



to the zero temperature case in Ref. . 

We treat the uniform matter and non-uniform matter consistently using the same RMF 
theory. Hence, all the resulting thermodynamical quantities are consistent and smooth in the 
whole range. Fig.6 shows the free energy per baryon F as a function of the baryon mass density 
Pb with various Yp at T = 0, 1, 5 and lOMeV. In the low density and finite temperature range. 
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the free energy is not significantly altered by variation in Yp, because the interaction between 
nucleons is very weak at low density. The inclusion of alpha-particles can reduce the free energy 
just before the formation of heavy nuclei where the alpha-particle fraction is some significant 
number as shown in Fig.4, so that the inclusion of alpha-particles can make the phase transition 
smoother. 

The pressure is calculated from F through the following thermodynamical relation, 



P 



nUdF/dns) . (39) 



We show in Fig. 7 the pressure as a function of pb with various 1^ at T = 0, 1, 5 and lOMeV. 
At low density and finite temperature, the pressure is mainly provided by the thermal pressure, 
which is proportional to the particle number density. When the alpha-particle fraction is a 
small number, the pressure is almost equal to n^T and independent of Yp. Just before the 
formation of heavy nuclei, the pressure curve shows some weak dependence on Yp due to the 
increase of Xa, then a sharp drop happens due to the formation of heavy nuclei. For high Yp, 
the pressure drops to the negative area. The pressure drop in the non-uniform matter phase 
has a strong Yp dependence. In high density region, all curves come back to the positive area, 
and have a rapid rise as the density increases. The behavior of the pressure at high density 
is determined predominantly by the contribution of the vector meson u and isovector-vector 
meson p as discussed in the Refs. |10|, |1^ . 



In Fig. 8, we show the results of the entropy per baryon as a function of pb with various Yp at 
T = 0.1, 1, 5 and lOMeV. The effects of the alpha-particles and heavy nuclei are clearly shown 
in the case of T = IMeV, the drops around 10^ g/cm^ are due to the increase of the alpha- 
particle fraction, while the drops around 10^^ g/cm^ are due to the formation of heavy nuclei. 
In the case of T = O.lMeV, the non-uniform matter phase appears beyond the lowest density 
limit, so the entropy has a strong Yp dependence due to the formation of heavy nuclei. On 
the other hand, the free gases outside heavy nuclei become dominant compositions of nuclear 
matter at high temperature, then the effect of the formation of heavy nuclei is weakened as the 
temperature increases. 

The properties of the EOS for neutron star matter were discussed in details in Ref. |p.0||. The 
neutron star masses and profiles using this EOS were shown and compared with other EOS in 
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the same reference. The properties of the EOS for hot dense matter in supernova core including 
leptons were discussed in Ref.[|13 . 



5 Conclusions and discussions 

The relativistic EOS of nuclear matter is designed for the use of supernova simulations in the 
wide temperature and density range with various proton fractions as listed in Table 2. We have 
adopted the RMF theory with the non-linear a and uj terms, which was demonstrated to be 
successful in describing the properties of both stable and unstable nuclei as well as reproducing 
the self-energies in the RBHF theory. We have worked out consistent calculations for uniform 
matter and non-uniform matter within the RMF theory. The Thomas-Fermi approximation is 
adopted to describe inhomogeneous nuclear matter, which can be modeled as a mixture of free 
neutrons, free protons, alpha-particles, and a single species of heavy nuclei. The inclusion of 
alpha-particles has significant effect around the formation of heavy nuclei and make the phase 
transition smoother. 

The relativity plays an essential role in describing the nuclear saturation and the nuclear 
structure it also brings some distinctive properties in the EOS comparing with the case 
in the non-relativistic framework 13]. Therefore, it is very interesting and important to 



study the astrophysical phenomena such as supernova explosion and the birth of neutron stars 
using the relativistic EOS through extensive comparison with the previous studies using the 
non-relativistic EOS. 
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Appendix A 

We list the definitions of tlie pliysical quantities tabulated in the EOS table 

• (0) temperature: T [MeV] 

The temperature T is written at the beginning of each block in the table, while we use 
mark 'cccccccc' to divide blocks with different T. 

• (1) logarithm of baryon mass density: logio(pB) [g/cm'^] 

• (2) baryon number density: ub [fm'^] 

The baryon number density is related to the baryon mass density as pb = rnuUB with 
niu = 931.49432Mel^ being the atomic mass unit taken from Ref. . 



(3) logarithm of proton fraction: \ogiQ(Yp) 

The proton fraction Yp of uniform matter is defined by 

~l~ 2?1q, Tip -\- ItIq, ( AfW 

Ub nn + np + Aria 

where Up is the proton number density, n„ is the neutron number density, na is the 
alpha-particle number density, and is the baryon number density. 

For non-uniform matter case, Yp is the average proton fraction defined by 

N 

yp = — (41) 

Nb ^ ^ 

where Np is the proton number per cell, Nb is the baryon number per cell; 

Np = I [np{r) + 2na{r)]d\ (42) 

J cell 



Nb = / [nn{r)+np{r)+Ana{r)]d\ (43) 

Jcell 

np{r) and n„(r) are the proton and neutron density distribution function given in Eq.(23), 
na{r) is the alpha-particle density distribution function given in Eq.(25). 

(4) proton fraction: Yp 
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(5) free energy per baryon: F [MeV] 

The free energy per baryon is defined as relative to the free nucleon mass M = 938.0MeV^ 
in TMl parameter set as 

f 

F=^-M. (44) 

(6) internal energy per baryon: Ei^t [MeV] 

The internal energy per baryon is defined as relative to the atomic mass unit as 

Eint ^ — - ruu. (45) 

(7) entropy per baryon: S [/cb] 

The entropy per baryon is related to the entropy density via 

S^^. (46) 
Ub 

(8) mass number of heavy nucleus: A 

The mass number of heavy nucleus is defined by 

^ = / [unir-) +np{r-)]ATir^dr, (47) 

J 

where n„(r) and ripir) are the density distribution functions in the Thomas-Fermi approx- 
imation as given in Eq.(23), Ra is the maximum between Rp and which is considered 
as the boundary of the heavy nucleus. 

(9) charge number of heavy nucleus: Z 

The charge number of heavy nucleus is defined by 

Z = / Tip ir) Anr^dr, (48) 
Jo 

(10) effective mass: M* [MeV] 

The effective mass is obtained in the RMF theory for uniform matter. In non-uniform 
matter case, the effective mass is a function of space due to inhomogeneity of nucleon 
distribution, so it is meaningless to list the effective mass for non-uniform matter. We 
replace the effective mass M* by the free nucleon mass M in the non-uniform matter 
phase. 
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(11) free neutron fraction: 

The free neutron fraction is given by 

X„ = (nri^°"*)/(nBV;e«) (49) 

where Vceii is the cell volume as defined in Eq.(22), = Vceii — is the volume 

outside heavy nucleus, n°"* is the free neutron number density outside heavy nucleus, 
while riB is the average baryon number density. 

(12) free proton fraction: Xp 

The free proton fraction is given by 

X^ = «V°"*)/(^BKeiO (50) 

with n""* is the free proton number density outside heavy nucleus. 

(13) alpha-particle fraction: X^ 

The alpha-particle fraction is defined by 

Xa = 47Va/(nsKeii) (51) 

where is the alpha-particle number per cell obtained by 

Na^ I n^{r)d\ (52) 

J cell 

na{r) is the alpha-particle density distribution function given in Eq.(25). 

(14) heavy nucleus fraction: Xa 

The heavy nucleus fraction is defined by 

Xa = A/iriBVceii) (53) 
where A is the mass number of heavy nucleus as defined in Eq.(47). 

(15) pressure: p [MeV / fm^] 

The pressure is calculated through the following thermodynamical relation 



nl{dF/dnB)\^^ (54) 
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• (16) chemical potential of neutron: /i„ [MeV^] 

The chemical potential of neutron relative to the free nucleon mass M is calculated 
through the thermodynamical relation 

^ln=[^{nBF)/^nn]T^^^ (55) 

here n„ = (1 — 1^) hb- 

• (17) chemical potential of proton: /ip [MeV] 

The chemical potential of proton relative to the free nucleon mass M is calculated through 
the thermodynamical relation 

^ip = [d{nBF)/dnp]^^^^ (56) 

here Up = YpUB- 

We have done the following check for the EOS table, 

(1) the consistency of the thermodynamical quantities, 

F + ^ = fir,{l-Yp)+fipYp, (57) 

Ub 

(2) the consistency of the fractions, 

Xn + Xp + X^+XA = l, (58) 

(3) the consistency of the relation between F, Eint, and S, 

F = Eint -TS + m^-M. (59) 

Gernerally these consistency relations can be satisfied within the accuracy of 0.001. Physical 
constants to convert units are taken from Ref . ||14|| . 
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Figure captions 

Fig.l. The distributions of neutrons, protons, and alpha-particles in the Wigner-Seitz cell at 
T = lOMeV, Yp = 0.3, and pb = are plotted by solid curve, dashed curve, and 

dot-dashed curve, respectively. In order to show the alpha-particle distribution clearly, we 
enlarge it by a factor 10. The parameters nf^, Ri and ti of the density distributions 
are defined in Eqs.(23) and (25), while Rceii represents the radius of spherical Wigner-Seitz 
cell. 

Fig. 2. The phase diagrams of the nuclear matter at 1^ = 0.01,0.1,0.2, and 0.4 in pB — T 
plane. Xa and Xa are alpha-particle fraction and heavy nucleus fraction, respectively. 
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The regions denoted by Xa > correspond to the non-uniform matter phase. The sohd 
curves show the boundary between the non-uniform matter phase and uniform matter 
phase, while dashed curves show the boundary where the alpha-particle fraction changes 
between X„ < 10"^ and > IQ-^ 

Fig.3. The phase diagrams of the nuclear matter at T = 0, 1, 5, and lOMeV in pB — Yp plane. 
The notations of X^, Xa and the boundary curves are the same as those in Fig. 2. 

Fig.4. The neutron fraction X„ (dotted curve), proton fraction Xp (dot-dashed curve), alpha- 
particle fraction Xa (dashed curve), and heavy nucleus fraction Xa (solid curve) as func- 
tions of the baryon mass density pb a-t Yp = 0.2 with T = 0, 1, 5, and lOMeV. 

Fig. 5a. The nuclear mass number ^4 as a function of the baryon mass density pb with the 
proton fraction Yp = 0.1, 0.2, 0.3, 0.4, and 0.5 at T = IMeV. 

Fig. 5b. The nuclear proton number Z as a function of the baryon mass density pb with the 
proton fraction Yp = 0.1, 0.2, 0.3, 0.4, and 0.5 at T = IMeV. 

Fig. 6. The free energy per baryon as a function of the baryon mass density pb with the proton 
fraction Yp = 0.01, 0.1, 0.2, 0.3, 0.4, and 0.5 at T = 0, 1, 5, and lOMeV. 

Fig. 7. The pressure as a function of the baryon mass density pb with the proton fraction 
Yp = 0.01, 0.1, 0.2, 0.3, 0.4, and 0.5 at T = 0, 1, 5, and lOMeV. 

Fig.8. The entropy per baryon as a function of the baryon mass density pb with the proton 
fraction Yp = 0.01,0.1,0.2,0.3,0.4, and 0.5 at T = 0.1 and IMeV, and with Yp = 0.01 

and 0.5 at T = 5 and lOMeV. 
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Table 1: The parameter set TMl for the lagrangian 



Parameter 


TMl 


M [MeV] 


938.0 


rua [MeV] 


511.19777 


[MeV] 


783.0 


rup [MeV] 


770.0 


9a 


10.02892 


9w 


12.61394 


9p 


4.63219 


92 [fm-^] 


-7.23247 


93 


0.61833 


C3 


71.30747 



Table 2: The ranges of the temperature T, proton fraction Yp, and baryon mass density pb 
the EOS table. The numbers of grid points in the ranges are also listed. Note that the cases 
T — and Yp — are covered with the listed ranges of the other parameters. 



Parameter 


Minimum 


Maximum 


Grids 


logio(T) [MeV] 


-1.00 


2.00 


31 




-2.00 


-0.25 


71 


logio(PB) [9/crn^] 


5.10 


15.40 


104 
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